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XIV. — The  Mathematical  Theory  of  Random  Migration  and 
Epidemic  Distribution.  By  John  Brownlee,  M.D.,  D.Sc. 

(MS.  received  May  16,  1910.  Eead  same  date.) 

The  general  theory  of  epidemic  disease  I have  already  considered  in  a 
communication  to  this  Society.*  In  that  communication  I showed  that 
the  course  of  epidemics  of  all  forms  of  infectious  disease  obeyed  certain 
very  definite  laws.  In  the  same  paper  it  was  also  shown  that  the  distri- 
bution of  epidemic  disease  in  a uniformly  populated  area  obeyed  a law 
essentially  similar.  Certain  reasons  were  given  why  the  normal  curve  of 

a:2 

error  y = y^e  might  be  expected  to  give  an  approximate  solution  in 

both  the  cases  considered,  but  why  the  distribution  actually  found  (type  iv.) 
should  be  the  common  form  was  not  at  all  clear.  I think,  however,  I have 
now  arrived  at  the  solution. 

The  distribution  of  an  epidemic  in  space  is  evidently  a problem  in 
chance.  If  there  is  an  infective  group  in  the  middle  of  a uniformly  dis- 
posed population,  then  the  distance  from  which  friends  come  to  visit  a sick 
person  or  the  distance  a sick  person  travels  while  developing  the  disease 
determines  the  subsequent  distribution  of  cases — a distribution,  therefore, 
obeying  some  law  on  the  average.  This  problem  has  since  been  attacked 
and  solved  by  Professor  Pearson  under  the  title  of  “ The  Problem  of 
Random  Migration.”  The  case  which  he  considers  refers  specially  to  the 
prevention  of  malaria,  which  is  now  known  to  be  spread  through  the 
agency  of  mosquitoes.  The  mathematical  theory,  which  is  very  complex, 
leads  to  the  determination  that  the  normal  surface  of  error  gives  a very 
close  representation  of  this  distribution.  For  epidemiological  purposes 
the  I’esult  is  quite  sufficiently  close.  To  make  the  matter  pei’fectly  clear, 
the  conditions  of  the  problem  solved  are  given  in  Professor  Pearson’s 
own  words ; — 

“ (1)  Breeding  grounds  and  food  supply  are  supposed  to  have  an 
average  uniform  distribution  over  the  district  under  considera- 
tion. There  is  to  he  no  special  following  of  river  beds  or  forest 
tracks. 

“(2)  The  species  scattering  from  a centre  is  supposed  to  distribute 
itself  uniformly  in  all  directions.  The  average  distance  through 
which  an  individual  of  the  species  moves  from  habitat  to  habitat 

* Proc.  Roy.  Soc.  Edin.,  June  1906. 
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will  be  spoken  of  as  a ‘ flight/  and  there  may  be  n such  ‘ flights  ’ 
from  locus  of  origin  to  breeding  ground,  or  again  from  breeding- 
ground  to  breeding  ground,  if  the  species  reproduces  more  than 
once.  A flight  is  to  be  distinguished  from  a ‘ flitter,’  a mere  to- 
and-fro  motion  associated  with  the  quest  for  food  or  mate  in  the 
neighbourhood  of  the  habitat. 

“(3)  Now,  taking  a centre,  reduced  in  the  idealised  system  to  a point, 
what  would  be  the  disti'ibution  after  random  flights  of  N indi- 
viduals departing  from  this  centre  ? This  is  the  first  problem. 
I will  call  it  the  ‘ Fundamental  Problem  of  Random  Migration.’ 

“ (4)  Supposing  the  first  problem  solved,  we  have  now  to  distribute 
such  points  over  an  area  bounded  by  any  contour,  and  mark  the 
distribution  on  both  sides  of  the  contour  after  any  number  of 
breeding  seasons.  The  shape  of  the  contour  and  the  number  of 
seasons  dealt  with  will  provide  a series  of  problems  which  may 
be  spoken  of  as  ‘ Secondary  Problems  of  Migration.’  ” 

The  proof  of  the  theory  given  by  Professor  Pearson  contains  also  im- 
plicitly the  proof  that  if  the  normal  surface  of  error  describes  the  distribution 
at  any  moment,  it  will  at  all  subsequent  times.  This  can  be  seen,  however, 

quite  easily  otherwise.  Thus  if  2/  = t>e  the  distribution  of  disease 

where  is  the  number  of  cases  per  unit  area  at  origin  and  Jx^  + is  the 
distance  of  any  point  from  the  origin,  then  the  amount  of  disease  at  any 

point  x',  y'  is  y^e  . This  element  gives  rise  to  a new  normal  surface 

_ (a:-g')2-Ky-y')2  ^ x"^+y"i 

y^e  2o-2  ^}^g  infection  at  x'y'.  If  these  several  surfaces 

be  integrated  from  — oo  to  -b  oo  with  respect  to  x and  y'  successively,  we 
get  the  new  distribution 

yoo  roo  -(x-x'^+^y-y'T-  g'2+^'2 

y = yA  I e 2<t2  yy.'^y' 

J ~'X)  J - CO 

a;2+2/2 

= y(pr<r‘^e  , 

SO  that  the  standard  deviation  o-  is  multiplied  by  s/2 — that  is,  that  the  slope 
of  the  surface  is  flattened.  In  other  words,  the  longer  the  disease  is  present 
in  a town  the  more  uniformly,  other  things  being  equal,  it  will  tend  to 
distribute  itself. 

At  the  present  point  it  will  perhaps  be  well  to  recapitulate  the  method 
by  which  it  was  considered  likely  that  the  normal  curve  of  error  should 
represent  the  course  of  an  epidemic  in  time.  If  there  be  an  amount  of 
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infectious  disease  a : if  such  element  infect  pa  persons,  and  if  q be  the  rate 
of  loss  of  infectivity  per  unit  time,  as  has  already  been  shown  to  he  probable, 


susceptible  persons  large)  is  a : ap  : apxp>q  : ap^qxpq^  : ap^q^xpq^,  etc., 

. . x-l 

or,  in  general,  if  x denote  the  unit  time,  y — aqo  q 2 , which,  as  q is 

necessarily  less  than  unity,  is  the  normal  curve  of  error.  If,  instead  of 
finite,  infinitesimal  differences  are  employed,  the  result  is  expressed  by  the 
formula 


That  the  normal  curve  is  to  be  taken  as  that  from  which  variation  is  to  be 
expected  both  when  the  space  and  time  distributions  of  epidemics  are 
examined  then  seems  clear,  and  it  remains  to  discover  in  what  manner  the 
natural  process  differs  from  that  so  far  developed  theoretically. 

It  is,  in  the  first  place,  to  be  noted  that  if  the  distribution  is  from  a 
central  area  instead  of  a point,  a disturbing  influence  on  the  shape  of  the 
curve  comes  into  play.  This  can  be  allowed  for  at  once.  For  purposes  of 
convenience  a two-dimensional  solution  is  given,  such  a solution  being  easily 
extended  to  three  dimensions  in  any  particular  case.  If  we  consider  the 
modification  of  the  normal  curve  produced  when  the  mosquitoes  start  from 
an  area  and  not  from  a point,  the  moments  of  the  resulting  distribution  will 
not  be  those  of  the  normal  curve,  that  is,  not 


the  number  of  each  group  of  infected  persons  (supposing  the  supply  of 


y = 


but 


or,  for  the  even  moments,  since  the  odd  moments  are  zero. 


so  that 


since  for  the  normal  curve  — 
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The  latter  is  always  less  than  3,  so  that  if  the  centre  of  distribution  is 
extended,  that  is,  if  the  original  mosquitoes  are  uniformly  dispersed  over  a 
space  bounded  by  two  parallel  lines,  the  subsequent  distribution  will 
resemble  a curve  of  type  ii.  rather  than  the  normal  curve.  But,  as  has  been 
said,  type  iv.  almost  uniformly  occurs.  Some  other  modification  is  therefore 
necessary.  This  may  be  found  in  the  fact  that  cr  is  not  constant.  In  the 
simplified  problem  of  Professor  Pearson  2cr“  = nV'  where  n is  the  number  of 
the  flights  and  I the  length  of  the  mean  flight.  The  curve  is  thus  given  by 
N 

„ ni2  j£  however,  vary,  we  do  not  get  the  normal  curve  at  all. 


y 


rnl^ 


but  a derivative.  The  frequency  of  any  value  of  f or  o-  may  be  taken  as 
given  by  /(<r)  when  the  limits  of  a are  a and  The  surface  of  distribution 
N 


derived  from  this  y = 


^ / 2<t2 


d(T.  F urther,  if  the  distribution  take  place 


from  a definite  area  as  above  described,  the  final  form  of  the  distribution  of 
the  organism  becomes 


fx+c 

' = / ' 
Jx-c 


cfixdx 


cr^ 


(1) 


when  (j)X  denotes  the  mode  of  distribution  in  the  area.  This,  if  the  forms  of 
(p{x)  and  {far)  are  known,  I take  to  be  the  fundamental  epidemic  or  random 
migration  equation. 

To  return  for  a moment  to  the  form  found  for  the  time  disti’ibution  of 
an  epidemic 


This  may  be  put  in  form 


y = ae 


y = ap^q'  • 

X\Ogp  + ix3  log? 


■ae  ^ log  7 '087 


(x+c)^_^_c 


= ae 


= ae 


changing;  the  origin  of  x 


(2) 


if  log  9 = - ^ and  = e. 

0--  log  q 


From  the  symmetry  of  the  epidemic  c must  in  general  be  a constant,  so  that 
we  have  the  relation 


or 


c log  q = log  p 


p = q<^ 


as  the  relationship  between  the  infectivity  and  the  rate  of  loss  of  infectivity. 
As  q is  less  than  unity  c will  be  negative  in  sign,  so  that  we  have  as  the 


266  Proceedings  of  the  Royal  Society  of  Edinburgh.  [Sess. 


general  rule  that  as  p increases  q decreases,  or  the  greater  the  infectivity 
the  more  rapidly  it  is  lost. 

In  the  case  of  the  epidemic  the  frequency  of  each  value  of  cr  is  quite 
2 

unknown.  As  log  q = and  as  q must  from  the  initial  assumption  lie 

between  0 and  l,the  limits  of  cr  are  evidently  — oo  and  0.  The  situation  of 
the  commonest  value  of  cr  is  therefore  between  these  limits,  and  must  in 
general  lie  nearer  zero  than  infinity.  In  choosing  an  arbitrary  form  for  cr 
so  as  to  get  an  approximation,  it  must  be  of  such  a form  that  (2)  will  be 
integrable  and  expressible  in  a form  suitable  to  calculation. 


r 

The  equation  obtained  by  integrating  (2)  is  y = I 


J __|2  ^ 

e fader.  It  is 


obvious  that  if  this  be  finite  when  a — 0 that  the  term  must  disappear, 


so  that  part  of  f(a)  must  be  e : the  other  part  may  be  taken  as  e 


*2' 


The  function  e 


has  a maximum  when  cr—  J ck.  The  constant  c is 


quite  unknown,  nor  does  it  seem  ascertainable  from  the  method  of  analysis, 
but  it  disappears,  and  as  k is  at  our  disposal  it  is  evident  that  the  maximum 
value  of  a can  be  placed  where  the  statistics  to  be  examined  demand.  We 
thus  have  as  a possible  form  of  the  epidemic  equation 


5 _^_0^ 

e ^dx. 


To  increase  the  variety  of  the  distribution  we  might  take  cr”e 


^2  0-2 


as 


representing  the  variation  of  cr : in  this  case  the  final  integral  becomes 

y = / o-"e  ^'^da. 

Jo 

The  epidemic  due  to  an  organism  instantsjneously  becoming  infective 
and  thereafter  losing  its  infectivity  at  a rate  corresponding  to  the  geo- 
metrical progression  should  at  least  approximately  fit  the  above  curve  if 
the  distribution  of  cr  has  been  at  all  closely  guessed.  When  we  turn  to 
Professor  Pearson’s  approximate  form  of  solution  of  the  random  migration 
problem  we  find  that  it  also  has  a term  with  a in  the  denominator.  The 
distribution  is  given  by 


N 


y = ¥r 


e 2(72 


where  N is  the  total  number  of  mosquitoes  starting  from  the  point  of  origin 
and  a = n being  equal  to  the  number  of  flights  and  I to  the  length  of 


1910-11.]  Mathematical  Theory  of  Kandom  Migration.  267 

the  average  flight.  Now,  if  a vary  from  0 to  x for  a first  approximation, 

_ CT 

its  frequency  might  be  represented  by  type  iii.  or,  say,  y = a-”e  *,  in  which 
case  the  migration  form  would  be 

/ crM’%  *• 

'2tt  Jo 

0-2 

This  form  is  very  intractable  as  a working  integral,  and  as  gives 


a somewhat  similar  series  of  curves,  the  latter  may  be  substituted  with 
much  simplification  of  the  calculation,  so  that  the  distribution  of  a species 
might  be  represented  by 

C I (T^"  + % ^'dor, 

Jo 

a form  easily  reducible  for  any  value  of  n.  This  is  the  same  surface 
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as  was  deduced  as  a possible  one  for  the  epidemic  time  wave.  What  has 
preceded  is  to  a certain  extent  hypothetical  in  so  far  as  the  assumption  of 
a form  for  the  frequency  of  o-  is  connected  ; otherwise  it  is  directly  deduced 
in  the  case  of  the  epidemic  from  the  consideration  of  a large  number  of 


typical  epidemics,  and  in  the  second  place  by  Professor  Pearson  by  a rigid 
mathematical  result  from  the  assumptions  already  referred  to. 

The  next  assumption  is  the  rate  at  which  infective  organisms  are  evolved. 
This  may  be  taken  as  represented  by  so  that  at  any  point  x the 

ordinate  will  be  represented  by  the  integral  of  the  curve  taken  between 
x — c and  x + c.  We  then,  finally,  obtain  as  the  form  of  the  epidemic 

fx+c  /•” 

y=yo\ 

Jo 

It  will  be  shown  later  that  good  agreement  of  theory  with  fact  is 
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obtained  i£  = 0 and  cf>{x)  = c (a  constant),  so  that,  in  the  first  instance,  the  form 


H — a 


x+c 


JO 

ilx  \ e 


will  be  considered. 

As  is  well  known 


3 X-  ^ 


the  positive  sign  referring  to  tlie  negative  branch  of  x and  the  negative 
sign  to  the  positive  branch. 

Hence  the  equation  of  the  curve  is 

_ 2c 

y = al\  — ake  ^ cosh  — from  x=Q  to  x — c, 
h 

- — , 2c 

ij  = ake  * sinh  --  from  x = c to  x = co  . 
k 


The  curve  is  symmetrical.  Examples  of  this  are  given  in  Diagrams  I. 
and  II.  When  these  are  compared  with  Diagrams  VII.  and  XXIII.  re- 
spectively of  my  former  paper,*  it  is  seen  how  much  better  the  fit  now 
obtained  is. 


Mathematical  Formulae  of  Curves  which  represent  possible 

Epidemic  Forms. 


Equations  and  formula  which  might  describe  more  or  less  approxi- 
mately epidemic  or  migration  forms  will  now  be  considered  in  detail.  The 
following  symbols  are  used  throughout 
(1)  Let  ?/ =/(x)  be  the  curve  of  distribution 

rcc 

A = I ydx,  area  of  the  curve 


where  the  axis  of  y is  taken  through  the  centre 
of  gravity 


c*cO 

x^ydx,  etc. ; 

— 00 


and 


xydx. 


Jo 

When  the  origin  is  not  in  a line  through  the  centre  of  gravity  the  moments 
are  denoted  by 


A*'  1!  M 2>  A*'  3> 

* Froc.  Roy.  Soc.  JSdin.,  vol.  xxvi.  pp.  491  and  507. 
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To  complete  the  symbols 


Q Q — 


/^4 


The  fundamental  forms  which  have  been  chosen  for  investigation  are 
/•=<= 

y = a\  a--'‘e  and  y - (if  o-"e 

Jn  Jo 

These  are  symmetrical  round  the  origin.  On  the  hypothesis  before 
mentioned  they  are  integrated  for  each  value  of  x from  x — c to  x + c, 
multiplying  each  term  by  a suitable  function  <p(x),  so  that  for  the  working 
equations  we  have, 


fx+c 

y = a (j>xdx  I o-^“e  ‘‘^do 

Jx-c  Jo 

fx+c  -^^-ya- 

t/  = a (jixdx  I a-"e 
jx-c  Jo 


(M) 


(N) 


(2)  The  first  of  these  forms  is  much  easier  to  evaluate,  and  in  addition 
gives  the  closest,  representation  of  the  facts.  Considering  it  in  the  first 
instance  and  assuming  that  (p{x)  = a (a  constant),  we  have  the  curve  of 
distribution  in  time  or  space  given  by 

fx+c  /-«> 

y = af  dxi  cr-”e  ■•^da-. 

' Jx-c  Jo 

This  form  implies  either  that  the  rate  at  which  infective  organisms  are 
given  off  is  constant  or  that  the  distribution  of  organisms  before  migration 
begins  is  uniform. 

o 

In  the  application  required  n is  in  general  zero.  The  areas  and 
moments  of  the  first  three  curves,  however,  are  given  in  the  following 
table ; — 


n = 0. 

71=  1. 

11=2. 

A 

ck’’  ijir 

ck“^  v/ 7T 

2ckP  sJtt 

Ic 

3k 

Ibk 

*'i 

2 

I 

16 

y-2 

k^  c2 
2 3 

/r2  + ^ 

3 ,,  c2 

k-i 

2 5 

Q ..4 

-A/r4+2c2/.;2  + ^ 
2 5 

sic*  + 6cV  + 

From  the  value  of  ^ obtained  from  the  statistics  can  be  easily  calcu- 

lated,  and  thence  k.  The  latter  value  may  be  compared  with  value  of  k 
deduced  from  the  value  of  i/^^,  and  the  nearest  value  of  n ascertained.  The 
equation  of  the  curve  corresponding  to  = 0 is  as  follows  : — 

- 2c  2x 

y = ak  — ake  cosh  * from  x = 0 to  x = c 

- 2*  2c 

Ic 


y = ake  * sink  — from  a;  = c to  z'  = co 
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Those  for  higher  values  of  n may  be  obtained  by  differentiating  the  above 
values  with  respect  to 

(3)  To  find  k and  c from  the  moments  is  easy,  but  for  convenience  a 
table  is  given  by  which  the  values  of  — may  be  obtained  when  the  value 

of  i.e.  — has  been  calculated. 

6 " ^ 

Table  showina;  the  values  of  - for  different  values  of 
^ c 


/.•■t  4/i-2  4 6/42*-^  6^-' 

9 


k 

1 

c 

63.  . 

Differences. 

■50 

■4314 

227 

•55 

■4541 

230 

■60 

■4771 

232 

■65 

■5003 

230 

■70 

■5233 

226 

■75 

■5459 

222 

■80 

■5679 

213 

■85 

•5892 

205 

■90 

■6097 

196 

■95 

■6293 

185 

TOO 

•6478 

346 

1-10 

■6826 

310 

V20 

■71.36 

275 

1^30 

■7411 

244 

V40 

■7655 

216 

T50 

■7871 

(4)  A special  case  arises  in  random  migration  if  the  area  to  be  invaded 
is  bounded  by  a straight  boundary  on  the  one  side  of  which  there  is  an 
infinite  field  uniformly  stocked  with  the  organism  which  is  migrating.  In 
the  case  of  animals  actually  migrating  this  form  has  been  found  by  Pro- 
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fessor  Pearson,  but  it  is  much  simpler  in  the  case  of  plant  life.  We  have  to 
deal  with  the  integral 

/•oo  — 2(a;4~^ 


y = a e 
Jo 

, - 2x 

Ka  -rr 

V = 

J i) 


di 


Application  of  this  will  be  given. 
(5)  Secondly,  the  form 


rx+c  /•«= 

y — a\  dx  \ o-«e  da- 

Jx-c  Jo 


may  be  considered.  In  this  the  distribution  of  frequency  of  a is  taken  as 
g-ji+2g-yo-_  ypg  moments  are  easily  calculated,  and  are  given  by 


y-n+2 

V — ^ ra  + 2 

^s/tt  7 

_c2  (n  + 3){n+-2) 

^2 


c^{n  + 3)(n  + 2)  3(re  + .5) {ji  + 4) («  + 3) (tH-  2) 


7“  7 

These  equations  are  solved  most  easily  if  c and  n be  eliminated.  This 
gives,  if  2 is  denoted  by 

9 .1  36.21  ^ 712.3^  6 A 1 , 7 9 2 , 6 .2  3 \ „ 

The  rest  of  the  solution  is  easy.  I have  not  been  able  to  arrive  at  the 
curve  without  mechanical  quadrature.  In  the  instance  in  which  it  has 
been  fitted  (Diagram  III.)  it  proves  exceedingly  unsuitable,  and  it  is  there- 
fore evident  that  the  form  discussed  cannot  represent  the  distribution  of 
a even  approximately. 


Asymmetry. 

(6)  Asymmetry  in  the  epidemic  curve  deserves  some  notice.  As  we 


found  before. 


logP 

log  q 


must  in  general  be  equal  to  a constant,  as  the  epidemic 


is  symmetrical.  If,  however,  jPMJhis  a function  of  a it  ma}^  on  account  of 

the  near  symmetry  of  the  epidemic,  be  assumed  that  it  can  be  expanded  in 
terms  of  a of  a rapidly  convergent  nature  so  that  all  the  terms  in  the 
expansion  except  that  of  a may  be  neglected. 
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The  expression 


becomes  on  this  hypothesis 


; I dx 

Jx-c  Jo 


In  this  case 


fx+c  /•=> 
a dx  I 
Jx-c  Jo 


- {X  - acr)2 


A =r  J irk^ 

fx  I — \J Irak 
, k'^ 

/^2  = 3+|-+a2F 

3 ~ 4 ^ + 2 ^ TTO^k^  + a.  J TTC^k 

k-'i  = + 12a2  + 3)  + ^V(l  + 2a^)  + 

D 

when  ^'2.  mV  Gtc.,  are  the  moments  round  the  origin,  which  gives 


M2  = o + 


k^ 


if  a2  = 0 


M3  = J \/nak^ 

3k'^  ,2  2 1 

M4  = -y  +^V  + -. 


For  the  form  of  the  curve 


the  solution  of  the  first  part  is  already  given.* 
The  second  is  equal  to 


H' ' {ix)dx* 
0 


This  method  for  accounting  for  symmetry  has  likewise  not  been  successful 
in  representing  the  statistics. 

(7)  When  the  distribution  is  symmetrical  around  a centre  it  is  evident 
that  we  have  to  deal  with  a similar  integration  of  Professor  Pearson’s  form 

_ ?• 

and  that  a distribution  of  y = ae  results. 


* A table  of  this  function  is  given  in  Jahnke  und  Emde,  Funktionentafeln,  p.  13.5. 


275 


1910-11.]  Mathematical  Theory  of  Random  Migration. 


This  can  be  fitted  by  radial  moments,  or,  in  general,  more  easily  by  first 
reducing  the  statistics  so  that  the  numbers  per  unit  zone  is  taken  as  a 
basis  of  calculation. 

The  distribution  may  also  be  summed  parallel  to  one  axis  and  compared 
with  a table  of  the  integral 

ale  ^ dx. 

Jo 

In  some  cases,  especially  when  an  epidemic  occurs  in  a locality  where  a 
disease  is  more  or  less  uniformly  endemic,  it  is  useful  to  have  the  means  of 
separating  the  epidemic  portion  of  the  disease  from  that  which  is  endemic. 
Such  a combination  is  very  common  in  the  case  of  such  diseases  as  scarlet 
fever,  enteric  fever,  diarrhoea,  etc.  Successive  approximations  are  necessary 
to  obtain  a solution,  but  a first  approximation  can  be  obtained  by  using  the 
normal  curve  to  represent  the  statistics  of  the  course  of  the  epidemic. 

Let  the  whole  amount  of  the  disease  be  as  represented  in  the  diagram 


base  and  the  epidemic  portion  by  the  curve  above.  In  general  the  parts 
of  the  curve  beyond  the  limits  may  be  neglected  as  only  affecting  the 
result  to  a small  extent.  Taking  the  middle  of  the  rectangle  as  origin, 
denoting  its  height  by  d and  its  length  21,  where  I is  known  and  d un- 
known, we  have  for  the  area  and  for  the  moments  of  the  rectangle  round 
an  axis  through  the  middle, 

A = 2W 


Afx„  = 

A/x^  = 


21^1 

2m 

5 ■ 


For  the  normal  curve  the  equation  is 


-{x-cT- 
V = ^0®  > 

and  its  area  and  moments  become 
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A'  = 2/o 

A/x'j  = y^sl  TTCrC 

_0-3  _ 

A/x  2~UoJ TT-^y-  + 2/o  TTO-C^ 

3 ~ //o  Vtt  -^  + //q  V • 


Adding  the  corresponding  portions  of  the  moments,  we  have  the  following 
equations  for  determining  the  four  unknowns,  d,  a,  c,  and  y^, 

A + Aj  = 2Zd  + y =Bj  . . . . . ■ {d) 

A/^\  = 2/o '/’TO'C  = Bg  ......  (b) 


A'  A ,-CT^  „ 

A //,  2 + A/Xg — + Z/o  v'tto-c^  = Bg 


. , i-<T^c  ,-crc° 

+yoAjTr-^- 


B, 


(«) 

(d) 


where  B^,  B2,  etc.,  denote  the  area  and  moments  obtained  from  the  statistics. 
By  using  (b)  the  equations  (a),  (c),  and  (d)  become 


From  a and  /3 


or 


which  gives  c. 

Hence  immediately  we  have 


?3  + 2W  = B, 
c 

p 2ld^  _ g 
— + n,C+  ^ Ug 

^B2  + ^B2  = B,  . 


0-2  B /2  72 

^B„  ~^  + B.,c  = B,  - — B,, 
2c  ^ 3c  ^ 3 


2B2C2  - (3Bg  - m^)c  + 3B^  - Bg/2  = 0, 


(a) 

(^) 

(r) 


(t2  = 2?4--c2 
Bg  3 

B, 

yo  = ~f= 

sItto-C 

2/cZ  = Bj-®2^ 

which  complete  the  solution. 

When  the  theory  just  described  is  applied  to  the  elucidation  of  epidemic 
processes  it  is  found  to  be  fairly  satisfactory.  It  gives  a very  good  fit  to 
many  epidemics,  as  will  be  shown.  It  is  also  capable  of  expressing  the 
moment  relationships  of  epidemics  satisfactorily.  Thus  it  is  generally 
found  that  /X4>3/X2:  this  is  explained,  but  the  form  also  allows  of  7x4  <87x9 
without  a change  of  theory. 
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Where  this  theory  fails,  however,  is  in  its  application  to  seasonable 
epidemics  such  as  enteric  fever  when  averaged  over  a number  of  years. 
This  might  be  expected.  Though  the  hypothesis  that  the  infectious 
organism  is  given  out  for  a definite  time  at  a constant  rate  may  be 
sufficiently  accurate  when  applied  to  a single  epidemic,  it  can  hardly  be 
expected  to  hold  when  an  average  of  years  is  considered.  An  early  or  a 
late  epidemic  of  the  disease  is  the  product  of  special  weather  conditions 
which  will  occur  very  seldom.  So  that  though  the  epidemics  occurring  at 
the  usual  time  will  tend  to  have  their  irregularities  smoothed,  the  beginning 
and  the  end  of  the  compound  curve  will  be  framed  on  a different  law  from 
the  middle.  An  endeavour  has  been  made  to  frame  an  approximate  formula 
to  meet  this,  but  none  has  been  found  fitting  the  facts. 

To  test  the  likelihood  of  a curve  fitting.  Professor  Pearson’s  method  is 
used  throughout.  To  do  this  the  actual  and  theoretical  numbers  are  differ- 
enced, the  difference  squared  and  divided  by  the  corresponding  theoretical 
number.  These  are  summed.  The  sum  is  denoted  by  The  table  for 
testino-  curve  fitting  * is  then  consulted,  and  the  value  of  P obtained. 

Thus  if  P = '8  it  signifies  that  in  testing  the  matter  over  again  a worse  fit 
might  be  expected  eight  times  out  of  ten,  so  the  fit  is  good.  The  opposite 
interpretation  is  given  when  P is  small,  say  = '2,  when  a worse  fit  might 
be  expected  only  in  two  out  of  ten  trials. 

Plague  in  Hong-Kono. 

The  epidemics  of  plague  in  Hong-Kong  are  typical  of  the  great  majority 
in  all  the  East.  They  are  nearly  bilaterally  symmetrical ; they  spring  from 


Table  showing  the  Number  op  Plague  Cases,  Actual  and  Theoretical,  in 
THE  Epidemics  op  Plague  in  Hong-Kong  in  1902  and  1904. 


1902. 

Number  of  Cases. 

1904. 

Number  of  Case.s. 

Actual. 

Theoretical. 

Actual. 

Actual  Figures  cor- 
rected so  as  to 
correspond  with 
the  Mean. 

Theoretical. 

March 

2 

3-8 

4 

3-6 

7-0 

A]U'il . 

27 

38-5 

40 

34-4 

35 

May  . 

157 

147-0 

135 

122-0 

126-4 

June  . 

194 

178-7 

194 

186-0 

180-5 

July  . 

131 

147-0 

96 

109-2 

126-4 

August 

50 

.38-5 

19 

29-3 

35-1 

September  . 

2 

3-8 

9 

11-3 

7-0 

Biometriha,  vol.  i.  p.  1.55. 
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a zero  line  where  the  disease  is  nearly  completely  absent.  The  only  point 
of  difficulty  in  treating  them  mathematically  is  in  settling  how  far  the 
calculation  of  the  moments  is  to  be  carried  into  the  interepidemic  period. 
In  general  good  fits  are  obtained  if  the  interepidemic  cases  are  neglected. 
Two  epidemics  are  shown  in  the  figures  given  in  the  table. 

The  fit  of  the  first  of  these  is  not  good  as  it  stands ; but  when  it  is  noted 
that  the  symmetry  with  regard  to  the  first  moment  is  produced  from  an 
unequal  distribution  with  respect  to  the  mean  of  a kind  that  the  sums  of 
the  equidistant  terms  on  both  sides  of  the  mean  are  nearly  equal,  we  may 
put  the  distribution  as  follows : — 


Actual. 

Theoretical. 

2 + 2 

4 

7-6 

27  + 50 

77 

77-o 

157  + 131 

288 

294-0 

194 

194 

178-7 

This  gives  — or  P = -28. 

Thus  the  variations  in  the  rise  cancel  those  of  the  fall.  This  method  is 
subject  to  criticism,  but  the  fit  of  any  individual  epidemic  can  hardly  be 
expected  to  be  good.  With  regard  to  the  epidemic  of  1904  the  fit  may  he 
said  to  be  good.  Here  = 5'96  or  P = '42. 

It  is  not  necessary  to  give  a large  number  of  examples  of  the  way  in 
which  the  present  theory  fits  the  facts.  In  many  instances  better  fits  are 
obtained  than  those  shown  in  my  previous  paper,  where  type  iv.  was 
found  to  closely  represent  the  epidemic  form.  Two  examples,  however, 
may  be  given,  that  for  the  smallpox  deaths  in  Warrington,  and  that  for  the 
milk  epidemic  of  scarlet  fever  in  Glasgow.  These  are  shown  in  Diagrams 
I.  and  ,11.  In  both  the  correspondence  of  the  facts  with  theory  is  very  close, 
much  closer  than  with  type  iv.  Before  leaving  this  part  of  the  subject  an 
example  of  the  use  of  the  distribution  of  a'^e~y'^  for  or  may  be  given.  This 
is  the  only  example  I have  thoroughly  worked  out,  but  in  a number  of 
others  I have  obtained  the  medium  value.  , In  none  is  there  any  evidence 
that  the  curve  obtained  has  any  resemblance  to  the  facts.  It  is  therefore 
very  improbable  that  can  represent  the  variation  of  a.  The  example 

illustrated  is  that  of  smallpox  in  London  in  1902.  As  will  be  seen  by 
referring  to  my  previous  paper,  type  iv.  gives  a quite  different  representa- 
tion from  the  curve  shown  in  this  case  (Diagram  III.). 

Examples  of  the  application  to  random  migration  will  now  be  given ; 
those  with  animal  forms  will  be  considered,  and  those  with  plant  forms 
thereafter.  With  the  former  it  is  difficult  to  secure  suitable  examples. 
Daphnia  pulex  was  used  in  many  experiments.  This  crustacean  does  not 
move  so  consistently  in  one  direction  as  others.  It  moves  by  jerks,  and. 
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moving  more  or  less  vertically,  fulfils  more  nearly  than  any  other  Professor 
Pearson’s  criterion  that  motion  in  any  one  direction  is  as  likely  as  in  any 
other.  It  has  also  the  advantages  of  being  not  specially  attracted  by  light 
and  of  being  more  or  less  opaque,  so  that  it  is  easily  photographed.  Cyclops, 
Littorina  rudis,  etc.,  were  also  used  as  is  described.  Of  plants  only  two 
species  were  investigated,  an  Oscillatorian  and  Aspargia  Uisp>ida. 

Daphnia  pulex. 

Experiments  were  made  in  two  ways  with  Daphnia  pulex.  In  some 
experiments  a large  number  of  the  crustaceans  were  placed  inside  a 
cylindrical  tube  in  a large  flat  dish  of  white  porcelain  and  then  liberated ; 
in  others  the  water  flea  was  allowed  to  take  up  its  position  as  it  liked  in 
the  dish.  It  usually  chose  to  distribute  itself  from  one  corner  along  one 
side  of  the  dish.  Examples  of  the  manner  in  which  this  happened  are 
given  below.  The  corner  being  an  impenetrable  boundary  may  be  taken  as 
representing  a centre  of  diffusion,  so  that  the  simple  fundamental  integral 
should  apply  and  the  grouping  should  conform  to  the  exponential.  This 
is  what  takes  place  (Diagram  IV.). 


Table  showing  the  Number  op  Daphnia  per  Unit  of  Length  along  the 
Margin  op  the  Plate  prom  one  Corner. 


Unit  of  Length. 

Actual  («). 

Theoretical. 

Actual  (6). 

Theoretical. 

0-1 

18 

18T 

23 

20-7 

1-2 

11 

lOT 

11 

13-6 

2-3 

3 

5-6 

7 

9T 

3-4 

4 

3-9 

6 

6-0 

4 — 5 

2 

2-5 

5 

3-9 

5-6 

1 

1-5 

1 

2-6 

(i-7 

1-1 

4 

1-7 

7-8 

1 

•7 

2 

1-1 

8-9 

1 

'75 

The  curve  is  given, 

X 

for  ((/)  by  ?/ = 18'e 


for  (6)byy  = 25’e 

In  the  first  example  the  fit  is  excellent  (y^  = 2'8  P = '9),  in  the  second 
not  so  good  (y^  = 6'34  P = -6).  The  want  of  fit,  however,  is  largely  due  to 
the  group  of  four  near  the  tail  end  of  the  distribution,  which  contributes  half 
of  the  divergence.  Except  for  this  group  a further  divergence  might  be 
expected  nine  times  in  each  ten  trials. 
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In  several  instances  the  Daphniee,  instead  of  grouping  themselves  along 
one  edge  of  the  dish,  arranged  themselves  as  if  one  corner  of  the  dish  were 
the  centre  of  attraction.  From  this  point  outwards  in  all  directions  of  the 
quadrant  their  numbers  diminished.  This,  however,  was  observed  only 
once,  when  the  light  was  good  enough  for  instantaneous  photography.  To 


/7 


20 

\ 



\ 

(6) 

\ 

s. 

\ 

/o 

N 

\ 

^ 

S 





count  the  numbers  of  the  organism,  concentric  circles  were  drawn  on 
the  plate  and  the  organisms  between  each  pair  of  concentric  circles 
in  the  quadrant  counted.  Some  difficulty  was  experienced,  however, 
at  the  greater  distances,  as  in  places  it  seemed  that  the  outlyers  of 
othei-  groupings  were  even  invading  the  periphery  of  this.  All  the 
organisms  were  counted,  however,  and  none  rejected.  The  figures  are  given 
in  the  following  table  : — 


i. 

Total 

Numbers. 

ii. 

iii. 

iv. 

Calculated 

Total. 

Number  per  Unit  Area. 

Actual. 

Calculated. 

Centre  qi^adrant 

19 

19 

18-4 

18-4 

1st 

,,  zone 

33 

11 

11-1 

33-3 

2nd 

37 

7-4 

6-78 

33-9 

3rd 

18 

2-67 

4-16 

29-12 

4th 

24 

2-67 

2-55 

22-95 

5th 

19 

1-72 

1-57 

16-27 

6th 

17 

1-31 

•96 

12'5 

7th 

11 

•73 

•59 

8-9 

178 

46-50 

46-11 

177-14 

It  is  easier  to  examine  this  grouping  when  the  numbers  are  reduced  to 
the  population  per  unit  area.  On  the  hypothesis  the  grouping  should  be 
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given  \>y  y = ae  * (par.  7 ) and  the  section  of  this  by  the  plane  of  y = 0 gives 

y = ae  k.  The  numbers  reduced  to  the  population  per  unit  area  are  given 

in  column  ii.  of  the  table.  If  we  take  the  area  and  first  moment  and 
integrate  from  cc  = 0 to  a;  = 8,  the  limit  to  which  a count  was  made,  we  get 

X 

y = 23’04e"2  056  ■ This  gives  or  P = -53,  a fit  as  good  as  could  be 

expected,  when  the  uncertainty  regarding  the  numbers  in  the  outer  zone  is 
remembered.  A large  part  of  the  value  of  y-  is,  as  before  noted,  due  to  one 
zone  alone.  In  all  these  groupings  there  seem  to  be  secondaxy  centres 
which  interfere  with  results  when  the  numbers  are  small. 

The  methods  of  dispersal  from  a centre  were  also  investigated.  F or  this 
many  Daphnim  (from  100  to  200)  were  placed  in  the  centime  of  the  dish 
with  a depth  of  water  of  about  of  an  inch.  They  were  contained  in  a 
cylindrical  tube  about  ^ inch  in  diameter.  When  the  level  of  the  water 
was  the  same  on  both  sides  of  the  tube  and  when  the  light  was  good  and 
the  camera  ready,  the  tube  was  removed,  the  dispersal  watched,  and  at  a 
suitable  moment  instantaneously  photographed.  The  photographs  of  course 
show  no  detail,  the  organisms  being  simply  marked  by  a paler  spot  on 
the  negative.  In  all  cases  a few  Daphnise  were  found  greatly  more 
energetic  than  the  rest.  These  were  generally  above  the  mean  size  and 
probably  represented  an  older  generation.  Tliey  were  so  exceptional  that 
they  possibly  should  be  rejected  from  the  statistics;  all,  however,  have 
been  included. 

The  experiments  took  much  time.  It  was  very  difficult  to  manipulate 
the  organism  without  damaging  a number  and  thus  introducing  a new 
factor. 

In  the  case  of  the  negative  from  which  the  following  table  is  made  the 
centre  of  the  group  was  found  by  counting  the  number  of  organisms  in 
each  half-inch  square  of  the  plate  and  calculating  the  mean.  This  being 
found,  circles  were  drawn  round  this  of  diameter  | inch,  1 inch,  1|  inches, 
etc.,  and  the  organisms  in  each  zone  counted. 

The  numbers  in  each  zone  are  given  in  the  table,  column  i. 

X 

The  curve  is  given  by  i/  = 39’2e'n369,  hence  y^  = 9'7,  which  gives  P = -4. 

In  this  case,  again,  one  zone  gives  a large  part  of  the  value  of  y^,  and 
it  is  also  to  be  noted  that  again  it  is  the  third  zone.  Another  fact  of 
interest  in  all  the  experiments  of  this  class  is  that  the  centre  of  the 
migration  is  not  the  original  centre  of  dispersal ; the  whole  mass  has 
moved  towards  the  light. 
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i. 

Total 

Numbers. 

ii. 

iii. 

iv. 

Total 

Numbers 

calculated. 

Number  per  Unit  Area. 

Actual. 

Calculated. 

Centre  circle  . 

28 

28 

27^8 

27^8 

1st  zone 

42 

14 

12^5 

375 

2nd  „ . . . 

28 

5-6 

5^8 

29^0 

3rd  „ . . . 

11 

1-57 

2^71 

19^0 

4th  „ . . . 

7 

•78 

L26 

1T3 

5th  . 

7 

■63 

■58 

63 

6th  . 

5 

•40 

■274 

3^6 

7th  „ . . . 

4 

■26 

■128 

U9 

8th  „ . . . 

2 

■12 

■06 

TO 

9th  „ . . . 

1 

■05 

■028 

■53 

The  dispersion  of  Daphnia  may  also  be  compared  with  the  formula 
- ^ 

y = ae  by  summing  parallel  to  one  axis.  A comparison  is  given  below 
when  the  Daphniee  were  so  counted.  The  numbers  are  as  follows : — 


Actual. 

Theoretical. 

2 

L9 

7 

9^4 

46 

39-5 

35 

39  5 

11 

9^4 

1 

L9 

giving  = 2'34  or  P = ’75. 


Cyclops. 

Cyclops  were  used  in  a few  experiments  but  found  very  unsatisfactory. 
Not  only  did  they  photograph  with  great^  difficulty  on  account  of  their 
transparency,  but  they  moved  on  dispersal  strongly  towards  the  light. 

A sectional  count  of  one  instance  is  given  (Diagram  V.).  It  is  not 
treated  as  a whole,  but  it  is  to  be  noticed  that  the  advancing  edge  is  very 
closely  given  by  the  exponential  curve  which  has  been  fitted  to  it. 

The  values  are  : — 

Actual.  Calculated. 

33  35-3 

20  19-2 

12  10-7 

7 5-9 

5 3-3 

2 1-8 


which  gives  y'^  = l’61  or  P = ‘9,  so  that  the  fit  is  very  close. 
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Othee  Crustacea. 

Some  copepods  were  experimented  with,  but  their  transparency  and 
rate  of  motion  unfitted  them  for  this  purpose ; they  were  also  too  strongly 
attracted  to  light.  In  addition  an  attempt  was  made  with  a fresh-water 


isopod,  but  a sufficiently  lai’ge  dish  could  not  be  obtained  to  make  the 
experiment  of  any  value. 

Littorina  rudis. 

Several  experiments  were  made  with  this  mollusc.  Numbers  were 
placed  in  various  vessels  and  their  movements  observed.  It  was  found 
best  to  put  them  at  the  bottom  of  a deep  glass  jar  with  a little  salt  water 


Distance  in 
inches. 

Numbers 

actual. 

Numbers 

theoretical. 

0-  -5 

25 

23-2 

■5-  1 

15 

16-3 

1 - 1-5 

13 

11-5 

1-5 -2-0 

6 

8-2 

2 0- 2-5 

2 

5-8 

2-5 -3-0 

5 

4-1 

3-0  - 3-5 

5 

3 0 

3-5 -4-0 

1 

2-2 

4-0 -4-5 

1 

1-4 

4-5 -5-0 

3 

ro 

5-0 -5-5 

1 

* i 

5-5 -6-0 

1 

'5 
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in  the  bottom  and  leave  them  to  climb  out ; this  they  did,  passing  the  edge 
and  falling  on  the  floor.  Latterly  a large  number  were  placed  in  the  jar 
one  day  and  next  day  those  on  the  side  of  the  jar  counted.  It  seems 
evident  that  we  are  dealing  with  the  tail  of  a moving  mass,  and  that  this 
tail  should  obey  the  law  of  the  exponential.  In  one  experiment  the  number 
of  Littorina  in  each  half  inch  of  the  jar,  measuring  from  the  top,  were  as 
shown  in  table  on  page  283. 

X 

This  gives  a distribution  y = 2'72e”2'W2’ 

As  it  stands  = 10'2,  so  that  P = '53,  but,  as  before  with  Daphnia,  40  per 
cent,  of  the  value  of  is  accounted  for  by  the  third  last  group.  If  this  be 
subtracted  P = ‘87.  In  either  case  the  lit  is  admissible. 
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Othek  Makine  Animals. 

With  regard  to  large  animals  few  experiments  can  be  recorded.  Some 
observations  were  made  on  flsh.  In  fish  migration  the  form  of  the  shoal 
should  be  capable  of  measurement.  In  this  case  the  tanks  in  the  Millport 
Biological  Station  seemed  to  aflbrd  some  chance  of  success.  One  of  the 
large  tanks  contained  a small  shoal  of  saithe  (Gadus  virens).  These  were 
photographed  several  times,  but  to  get  sufficient  light  for  an  instantaneous 
plate  was  difficult.  In  one  such  photograph  the  shoal  is  on  the  point  of 
turning.  The  symmetry  is  remarkable,  the  numbers  from  the  left  to  the 
right  in  each  unit  of  length  being  as  follows : — 

1,  2,  2-5,  6,  6,  2-5,  1,  1.  Total  22. 

The  numbers  are  too  small  to  allow  of  differentiation  of  the  type  of  curve, 
as  they  can  be  fitted  either  to  the  exponential  curve  described  or  to  the 
normal  curve  of  error. 
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Data  for  the  other  fish  are  lacking  with  the  exception  of  herring.  The 
weekly  takes  at  given  stations  afibrd  some  guide,  but  not  much.  One 
selected  from  the  Fishery  Board  Report  of  the  take  of  herring  ofi'  Ballantrae 
is  illustrated  in  the  diagram  (Diagram  VII.).  Here  the  elevation  at  the 
beginning  of  the  ascent  prevents  accurate  curve  fitting,  but  the  general 


correspondence  is  good.  It  was  hoped  that  with  the  figures  at  a fixed 
point  some  better  results  might  be  obtained,  and  the  Fishery  Board  very 
kindly  furnished  me  with  the  figures  of  the  weekly  takes  in  the  trammel 
nets  at  Ballantrae.  But  these  figures  proved  insufficient  to  determine  the 
shoal  form  ; in  no  case  wei’e  the  curves  continuous.  Factors  such  as  the 
end  of  the  open  time,  storms,  etc.,  interfered  with  the  returns  to  such  an 
extent  that  it  was  not  possible  to  use  the  figures.* 

* Herring  catch,  Ballantrae,  1903,  May-August,  Twenty-fifth  Ann.  Rep.  Fish.  Board  Sc., 
p.  174. 
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OSCILLATORIA. 

A species  of  Oscillatorian  (not  identified)  was  also  experimented  with. 
This,  when  put  in  a mass  in  the  bottom  of  a test  tube,  started  to  climb  the 
wall  to  the  surface.  After  the  migration  had  taken  place  so  that  the 
beginning  and  end  of  the  mass  was  fairly  clear,  the  test  tube  was  carefully 


emptied  and  filled  with  melted  gelatine.  When  this  was  set  the  number  of 
filaments  per  millimetre  of  length  on  the  test  tube  was  counted  under  the 
microscope.  The  numbers  in  each  space  of  four  millimetres  were  as 
follows : — 

2,  20,  42,  80,  125,  190,  168,  154,  164,  128,  70,  28,  10,  2. 

When  fitted  to  a migration  curve  the  correspondence  is  as  illustrated  in  the 
accompanying  diagram  (Diagram  VIII.).  The  fit  is  not  good,  but  the  shape 
of  the  curve  suggests  that  it  is  made  up  of  two  and  not  one  migration 
system. 
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Apargia  hispida. 

In  the  Belvidere  Hospital  there  are  long  straight  terraces  adjacent 
to  the  Clyde  constructed  as  in  the  diagram  (Diagram  IX.).  It  is  the 
custom  to  mow  the  level  parts  of  the  terraces  with  a lawn  mower. 
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No  plant,  therefore,  of  any  high  habit  can  ever  seed  there.  It  can 
only  grow  herbaceously.  On  the  bank  between  the  terraces,  however, 
such  plants  can  develop,  and  a considerable  stretch  of  the  bank  is 
thickly  overgrown  with  Apargia.  The  seeds  of  this  scatter  over  the 
lower  and  higher  terraces,  when  they  germinate  and  form  plants  in 
the  grass. 

This  case  can  easily  be  considered  as  one  where  there  is  an  infinite 
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uniform  source  on  one  side  of  a straight  line,  and  therefore  may  be 

X 

represented  by  y — ae^k' 

To  obtain  the  form  of  distribution  the  ground  on  both  terraces  was 
lined  with  parallel  lengths  of  string  two  feet  apart  for  a distance  of  fifteen 
yards,  and  the  number  of  plants  in  each  rectangular  space  counted.  The 
numbers  are  given  in  this  table ; — 


Number  oi’  Plants. 


Distance. 

Upper  Terrace. 

Lower  Terrace. 

0-  2 feet 

140 

201 

2-  4 „ 

47 

157 

4-  8 „ 

39 

99 

8-10  „ 

24 

50 

10-19  „ 

17 

22 

12-14 

10 

20 

14-16  „ 

9 

16 

16-18  „ 

5 

4 

18-20  „ 

7 

12 

20-22  „ 

3 

6 

22-24  „ 

4 

3 

24-26  „ 

3 

3 

26-28  „ 

1 

3 

28-30  „ 

3 

4 

30-32  „ 

1 

32-34  „ 

1 

Total 

312 

605 

The  equation  of  the  theoretical  curves  are : — 

X 

Higher  terrace  y = l‘l'le~ 2-45S' 

X 

Lower  „ 1/ = 248e”2-445- 

The  two  areas  show  a practically  identical  distribution,  with  the  excep- 
tion that  the  seed  spreads  to  twice  the  extent  on  the  lower  than  on  the 
higher  terrace.  The  nature  of  the  fit  is  shown  on  the  diagram,  and  as  the 
soil  of  the  locality  is  all  forced,  clay  coming  to  the  surface  in  patches,  and 
drainage  being  very  irregular,  it  is  as  good  as  might  be  expected. 


Conclusions. 

(1)  The  general  principles  which  underlie  both  epidemic  distribution  in 
space  and  time  and  random  migration  are  identical. 

(2)  Both  can  be  deduced  almost  directly  from  the  laws  of  chance  through 
assumptions  which  have  considerable  a priori  probability. 
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(3)  It  is  found  in  both  cases  that  the  exponential  curve  can  be  taken  as 
giving  a fair  representation  of  most  of  the  facts.  The  manner  in  which 
this  curve,  however,  appears  is  somewhat  different  in  the  two  instances. 
It  is  doubtful  why  this  formula  should  express  both  results. 

Note. — The  experimental  work  recorded  in  this  paper  was  done  chiefly 
at  the  laboratory  of  the  Millport  Marine  Biological  Station,  and  I desire 
to  thank  the  Superintendent  for  much  assistance. 


{Issued  separately  February  25,  1911.) 


